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Abstract 

Convergence of the solutions of nonhomogeneous linear singularly perturbed systems to that 
of the corresponding reduced singular system on the half-line [0, oo) is considered. To include 
the situation on a neighborhood of initial instant, a boundary layer, a distributional approach 
to convergence is adopted. An explicit analytical expression for the limit as a distribution is 
proved. 
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1 Introduction 

A rational motivation to study singular linear system, 

Ex = Ax + Bu (1.1) 

with singular matrix E, is that it is an evident simplification of the singularly perturbed systems 

E{e)x = Ax + Bu (1.2) 

for a "small" parameter e (may be of vector form), where E{e) is nonsingular and tends to E as 
e ^ 0. The system (1.2) arises naturally from, for example, coupling subsystems with "slowly" and 
*This work was supported by the Harbin Institute of Technology Science Foundation under grant HITC200712. 
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"fastly" varying states respectively, optimal linear-quadratic regulator with cheap control, etc. For 
detail, see [l]-[8]. For a specific system analysis or synthesis problem, the effectiveness of the above 
simplification relys on "approximate extent" between the solution to the problem for (1.1) and that 
for (1.2). Partially for characterizing "approximate extent" in the singular perturbation analysis, 
some interesting topologies are introduced. See, for example, [8], [9] and the references therein. 
In this paper, we are interested in the following singularly perturbed initial value problem 

N{e)x{t) = x{t) + f{t), t>Q 

[i.6) 

x{0) = xo, 

and the corresponding reduced one 

Nxit) = xit) + fit), t>0 ^^^^ 

X{0) = Xq, 

Here A'"(e) e R"^" is nonsingular for e 7^ and tends to N, a nilpotent matrix, as e ^ 0. The index 
of nilpotency of N is denoted by q, i.e., 

q^mm{k : k>l,N'' ^0}. (1.5) 

The nonhomogeneous term / is a g — 1 times continuously differentiable function mapping = 
[0,+oo) to R". Under a regularity assumption, the singular system (1.1) can be transformed into 
two subsystems through Weierstrass decomposition [10]. One has the form of the normal hnear 
system which has trivial relationship to the corresponding perturbed ones, and another is of the 
form (1.4). For more detail of background, see [2]. For general initial conditions ("inconsistent initial 
conditions"), the problem (1.4) has no solution in the sense of classical differentiable function, and 
the corresponding physical system exhibits impulsive behavior [1]. Thus some generalized solutions 
are adopted for the problem (1.4). Recently [11]-[13], an exphcit distributional solution of (1.4), 

<t) = - Eto N'f^Ht) - W=\ S^'-'\t)N'^ {xo + EVo iV7«(0)} ,t>0, (1.6) 

is obtained by Laplace transform. So in what sense and whether the solution of (1.3) given by 

x,{t) = eM{N{e))-H}xo + fl exp{(7V(e))-^(i - T)}{N{e))-'f{T)dT, t > 0, (1.7) 

a classical function mapping M+ to M"', can be approximated by the distribution (1.6) becomes 
interesting. 
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Works [6] and [7] have the same concern, but they only considered natural response (i.e., the 
solution for / — 0). The forced response (i.e., the solution for xq = 0) of (1.4) also contains 
impulse term at initial instant according to (1.6). So the convergence in a neighborhood of t = 0, 
a "boundary layer" (region of nonuniform convergence, see [6], [14]), for the forced response also 
appeal to a distributional approach. This motivates a generalization to the results in [7] to include 
the nonhomogeneous case. For other related works, see [3], [4], [14], [15] and the references therein. 

2 Notations and Definitions 

We review some notations and definitions in distribution theory [16]. Let Cg^(M, M") be the space 
of infinitely differentiable functions from R to with compact support. There is a topology on 
it [16], and then the distribution space is defined as the dual space C^(R, R")'. So a distribution 
w e Cg^(]R, M")' is a linear continuous functional on C^(]R, R"). The value, a real number, of w on 
A e C^(M,M") will be denoted by {w, A) . The Dirac delta distribution 5 E C^(M,M)' is defined by 
{5, A) = A(0) for VA e Cg'(R,R). For any distribution w e Cg'(R,R)', its k-th order distributional 
derivative V e C^(R, R)' is defined by 

(Pf ^, a) = (-1)^ A^'^)) (2.1) 

for VA e Cg'(R,R), where A^''^ denotes the k-th order usual derivative. Let >Cioc(R,R'') denote the 
set of all locally Lebesque integrable functions from R to R". The embedding map S : £ioc(R, R") 
^ C:^(R,R")' is defined by {Sz, A) = ^(t)TA(t)dt for G £ioc(R,R") and VA G C^(R,R"). 
Here z{t)^ represents the transpose of z{t), and the integral is in the sense of Lebesque. We do not 
distinguish z and £z in following. Lastly, let C'^(R+,R") denote the set of all A;-times continuously 
differentiable functions from R+ to R", which can be seen as a subset of >Cioc(R, R") naturally. 
Now we cite the definition of convergence of distribution sequence [16]. 

Definition 2.1 Given sequence {zi}°Zi C Cg^(R,R")' and z e Cg'(R,R")', then {zi}°Zi is said to 
converge to z in C^(R, R")', denoted by limj^oo^i = if for every A e C^(R, R"), 

limj^oo {zi, A) = (z. A) . 

For convenience and without loss of generality, we consider discrete perturbations 

Nix{t) = x{t) + f{t), t>0 ^^^^ 
x{0) = xo, 
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where Ni is nonsingular, and 

linii^ociVi = (2.3) 

The solution of (2.2) is 

Xi{t) = exp{Ar.-it}xo + /* exp{iVr'(t - r)}7VrV(^)c?^, t > 0. (2.4) 

Then we need to explore, in the sense of Definition 2.1, the convergence of the solution sequence 
to the solution (1.6). 

In following, except 5^''\ the k-th order derivative notation z^'^'' will always be in the ordinary sense 
according to pointwise differentiation. In the case z e C'^(R+,R"), notation z'^'^^O) is understood as 
that from right hand. We always assume / G C''~^(R+,M") in this paper, where q is the nilpotency 
index of N, to guarantee the distributional solution having the expression (1.6). 



3 Uniqueness 

For a perturbation manner given by {A^i}^i, the solution sequence {xi{t)}'^i may not converge. 
But we will prove that if it does, then the limit must be the solution (1.6) of the reduced system 
(1.4), not dependent of the perturbation manner. This generalizes Theorem 2 in [7]. 

Lemma 3.1 [16, p. 21]Let z e C*=(M+,R'^). Then we have 

Vfz = z^'^ + E -lo S^h^''-^-'\0). (3.1) 

Note that, according to the convention in Section 2, the precise meaning of (3.1) is 

£:(z) = £{z^''^) + s^^z^''-'-^\o). 

Lemma 3.2 Xi{t) e C«(R+, W) and /or m = 1, 2, . . . , 

xt^\t) ^ Nr^Xiit) + El^i Nr^f^^-^tl t > 0. (3.2) 
Proof. Firstly, we prove the case m = 1. 

= N-'x,{t) + N-'m. (3.3) 
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Secondly, supposing that the case m holds, we prove the case m + 1. Differentiating two sides of 
(3.2) gives 

= Nrxf^it) + Elii N.^f^-^'-\t). (3.4) 

Substituting (3.3) in (3.4) gives the result immediately. ■ 
Combining Lemmas 1 and 2, we have 

Lemma 3.3 For k — 1,2, . . . ,q, we have 

Vfx, = Nr^x, + Ef=i 7Vr7^'-'^ + E K 8^'^ (n7^'-'-'^x, + Et;"' A^r'/('-'-^-'^(0)) ■ (3.5) 

Lemma 3.4 /i^, p. 28] Let {zi}^^ C C^(M,M'*)' and^ e C^(M,M")'. IfYim^^^Zi = z, then for every 
k>l, 

lim Vfzi = Vfz. 

i— »oo 

Theorem 3.1 //{a^i}^! converges, then limj_»oo2^i = 2;. 

Proof. Let k = q, the index of nilpotency of N, in (3.5). Multiplying two sides from left by Nf gives 

- es ^^i/^^^ + Es; s^^^Nr {xo + Er=v ^vr/^-^^ (o)) . (3.6) 

Letting i ^ oo and noting that iVj — > A^, we obtain 

7V« lim V^^^x, = lim + ES A^V^'^ + EfJ (^o + E'Jo' ^"/^™^ (0)) 

from Lemma 3.4. Noting that ^ and N^+^ Hm=q-2~j+i = 0' ^e have 

hm X, ^ - N'f^ - E?:o ^^'^N^^' (^0 + Erio ^"^/^"^^ (0)) 

I— »oo •' 

= - EL^J N^f^'^ - ml s^'-'^N' (^0 + Erio ^"V^'"^ (0)) . 

This completes the proof. ■ 

4 Convergence 

In this section, we will give a condition on perturbation to guarantee convergence. An example 
satisfying the condition shows the existence of convergent perturbation. This gives a generalization 
to Theorem 1 in [7]. 



Lemma 4.1 // the number sequence 

{/o II^M* ^\\dt, 1,2,...} is bounded for some k > 0, 

and f e C'i+''{R+,W) n £\R+,W) then {xi}^^ converges. 

Proof. Under the boundedness assumption, the sequence {N?'^^'^'' e^i ^*xo}^i converges to in the 
sense of Definition 2.1 by Lemma 1 in [7]. Let h e C^(R,R") with \\h{t)\\ < C for e R. Since 

[Nr^'iy-'^'-'^N-^Mdr, h] 

< it^ mm ■ (/o ■ ii/wim-) dt 

< ir 11/(^)11 • {lTmr'e''^'^'-^^\\ ■ Whmdt) dr 

< ir\\fir)\\dr-\\N!\\CJt-\\N^e'^^''\\dt 


by the assumptions (note that ||iVf|| = 0), the sequence {A^f+^+*^ /* e^r'(*-^)Arri/(r)rfr : 

i = 1,2, . . .} converges to in C^{R, M")' also. So we have 

converges to in C^(]R, M")'. By Lemma 3.4, we have 

hm Pj+'+'=)(7V/+^+'=Xi) = hm iVf+^+'^Pi^+'+'^Xi = 0. (4.1) 

i— >oo i— >oo 

On the other hand, since / e C«+'=(R+, R"), we have 

= x, + Y:to'^'^~'Nif'^ (4.2) 



like (3.6). Prom (4.1) and (4.2) we see the existence of liuii^^Xi and 

hm X, = - EK'^'^"' N^f'^ - E?-o'^'^"' S^^^N^^' (xo + Yt-l^'^-'-' 7V-/M(o)) . 

i— >oo \ / 

Noting that A?"^ = 0, we see that it equals x by (1.6). ■ 

We intend to weaken the higher differentiability requirement for / e C^''"'^(R+, R") in Lemma 4.1. 
Again, we note that / is always assumed in C^~^(R+,R"). 

Lemma 4.2 Suppose f e >C^(R+,R'^). // the number sequence {Jq°° \ \N^e^i ^*\\dt, i — 1,2, ...} is 
bounded for some k >0, then {xi}'^^ converges. 
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Proof. We only prove the result in the case A; = 0. That for A; > 1 can be proved by some slight 
modification. Note that / e C9-i(R+,M") but maybe / ^ C«+°(R+,M") = C5(R+,M'*). 
Diff'erentiating two sides of (3.6) gives 

mV^r'\, = V^x, + NlV^f^'^ + Eg S^^'^Nf (x, + Er=V A^r/^-^Ho)) . (4.3) 

Noting that Xi e C^(K+, M") and e C^(M+, R") for Z = 0, 1, . . . , ? - 2, it follows from Lemma 3.1 
that 

V^Xi = Xi + 6-XQ 

and 

P,/(0 = + 5-/(0(0) 
for / = 0, 1, . . . , g — 2. Substituting in (4.3) gives 

+ ES A^'^^ (xo + Er=V A^r/^"^^ (0)) 

Then we have 

+ ^.-^^50+1)^+1+1 (a:o + Er=VA^r/^'"nO)) - A^f/^^-'Ho). 

Noting that 

lim Ar/Pd/(<?-i) = mV^f^'i-^^ = 0, 

i— >oo 

the remainder thing is similar to the proof of Lemma 4.1. ■ 

We need to weaken the integrability requirement / e >C^(R+,R'^). 

Lemma 4.3 For any 6 > 0, there exists h G C«-i(R+,R") n£i(R+,R") such that 

hit) ^ fit), yt<h. (4.4) 
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Proof. One can construct a (unique) polynomial P{t) of degree {2q — 1) such that 

pW(6) = P«(6 + l) = 

for A; = 0, 1, . . . , g — 1 (see [17, p. 88]). Then we define 



fit), ifO<t<b, 
P{t), iib<t<b + l, 
0, ifi>6 + l, 



which satisfies the requirement. 



Theorem 4.1 If the number sequence 

Uo I l^f ^\\dt, 1,2,...} is bounded for some k > 0, 

then {xi}'^i converges. 

Proof. Arbitrarily choose h e C^(R,R"). Then we have 

h(t) = 0, > 6 

for some b > 0. Let fb G C«-^(M+,M") n £^(M+,R") with (4.4). Then by Lemma 4.2, the sequence 

y,{t) = e^^'^xo + /oe^"(*-^)A^rV6(r)dr, t > 

converges to 

y{t) - - EJ=o ^Vi'^(t) - nil s^'-'Hm'^ {x-o + eVo ^Vf (0)} ,t>o 

in the sense of Definition 2.1. By direct computation we can get 

{xi, h) = {yi, h) , i = 1, 2, . . . 

and 

{x, h) = (y, /i) . 

Therefore limj^oo {^i, h) = {x, h) , and this completes the proof. ■ 

Example 4.1 Set Ni = N - \l , i = 1,2, .. .. Then {/+~ | jA/'f e^^'*| i = 1,2, . . .} is bounded for 
some /c > (see Lemma 2 in [7]). So according to this perturbation manner. Theorem 4.1 guarantees 
that the solution sequence {xi\'^i of the perturbed systems (2.2) converges to the solution x of the 
singular system (1.4). 



5 Conclusions 



As an idealized model, the nonhomogeneous singular system can approximate some singularly 
perturbed systems well in a sense of distribution theory. A future work is to give some condition 
easy to verify on perturbations to guarantee convergence. 

References 

[1] G. C. Verghese, B. C. Levy, T. Kailath, "A generahzed state-space for singular systems," IEEE 
Transactions on Automatic Control, Vol. AC-26, No. 4, pp. 811-831, 1981. 

[2] L. Dai, Singular Control Systems, Springer- Verlag Berlin, Heidelberg, 1989. 

[3] P. V. Kokotovic, R. E. O'Malley, Jr., and P. Sannuti, "Singular perturbation and order reduction 
in control theory — an overview," Automatica, Vol. 12, pp. 123-132, 1976. 

[4] P. V. Kokotovic, J. J. AUemong, J. R. Winkelman, and J. H. Chow, "Singular perturbation and 
iterative separation of time scales," Automatica, Vol. 16, pp. 23-33, 1980. 

[5] B. A. Francis, "The optimal linear-quadratic time- invariant regulator with cheap control," IEEE 
Transactions on Automatic Control, Vol. AC-24, pp. 616-621, 1979. 

[6] B. A. Francis, "Convergence in the boundary layer for singularly perturbed equations," Auto- 
matica, Vol. 18, No. 2, pp. 57-62, 1982. 

[7] D. Cobb, "On the solutions of linear differential equations with singular coefficients," Journal 
of Differential Equations, Vol. 46, pp. 310-323, 1982. 

[8] D. Cobb and J. Eapen, "High-gain state feedback analysis based on singular system theory," 
SIAM J. Control Optim. Vol. 44, No. 6, pp. 2210-2232, 2006. 

[9] D. Hinrichsen and J. O'Halloran, "Limits of generalized state space systems under proportional 
and derivative feedback," Math. Control Signals Systems, Vol. 10, pp. 97-124, 1997. 

[10] F. R. Gantmacher, The Theory of Matrices, Vol. 2, Higher Education Press, Beijing, 1955 
(Chinese translation from Russian). 



9 



[11] p. C. Miiller, "Remark on the solution of linear time-invariant descriptor systems," PAMM Proc. 
Appl. Math. Mech., 5, pp. 175-176, 2005. 

[12] Z. Yan and G. Duan, "Time domain solution to descriptor variable systems," IEEE Transactions 
on Automatic Control, Vol. 50, No. 11, pp. 1796-1799, 2005. 

[13] Z. Yan, "Geometric analysis of impulse controllability for descriptor system," Systems and Con- 
trol Letters, Vol. 56, No. 1, pp. 1-6, 2007. 

[14] R. E. O'Malley, Jr., "A singular singularly-perturbed linear boundary value problem," SI AM J. 
Math. Anal. Vol. 10, No. 4, pp. 695-708, 1979. 

[15] S. L. Campbell, Singular Systems of Differential Equations, Pitman, London, 1980. 

[16] I. M. Gelfand and G. E. Shilov, Generalized Functions, Vol. I, Science Press, Beijing, 1965 
(Chinese translation from Russian). 

[17] X. He, Y. Su, and X. Bao, A Short Course in Numerical Mathematics, Higher Education Press, 
Beijing, 1980 (in Chinese). 



10 



